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(57) A sampled amplitude read channel for mag- 
netic disk recording which asynchronously samples (24) 
the analog read signal, adaptively equalizes (B103) the 
resulting discrete time sample values according to a tar- 
get partial response, extracts synchronous sample val- 
ues through interpolated timing recovery (BIGG), and 
detects digital data from the synchronous sample val- 
ues using a Viterbi sequence detector (34) is disclosed. 
With interpolated timing recovery, the equalizer filter 
and its associated latency are removed from the timing 
recovery loop, thereby allowing a higher order discrete 
time filter and a lower order analog filter. 
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Description 

FIELD OF INVENTION 

5 The present invention relates to the control of magnetic disk storage systems for digital conoputers. particularly to 
a sampled anrplitude read channel that employs asynchronous sarrpling of an analog read signal, adaptive discrete 
time equalization, and interpolated timing recovery 

GROSS REFERENCE TO RELATED APPLICATIONS AND PATENTS 

10 

This application is related to other co-pending U.S. patent applications, namely application serial numbers 
08/440,515 entitled "Sampled Amplitude Read Channel For Reading User Data and Embedded Servo Data From a 
Magnetic Medium." 08/341,251 entitled "Sampled Amplitude Read Channel Comprising Sample Estimation Equaliza- 
tion, Defect Scanning, Channel Quality, Digital Servo Demodulation, PID RIter for Timing Recovery, and DC Offset 

15 Control," 08/313,491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording," and 
08/533,797 entitled "Improved Fault Tolerant Sync Mark Detector For Sampled Amplitude Magnetic Recording." This 
application is also related to several U.S. patents, namely U.S. Pat. No. 5,359,631 entitled "Timing Recovery Circuit for 
Synchronous Waveform Sampling," 5,291,499 entitled "Method and Apparatus for Reduced-Complexity Viterbi-Type 
Sequence Detectors," 5,297,184 entitled "Gain Control Circuit for Synchronous waveform Sampling," 5,329,554 enti- 

20 tied "Digital Pulse Detector." and 5,424,881 entitled "Synchronous Read Channel." All of the above-named patent appli- 
cations and patents are assigned to the same entity, and all are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

25 In magnetic storage systems for computers, digital data serves to modulate the current in a read/write head coil In 
order to write a sequence of corresponding magnetic flux transitions onto the surface of a magnetic medium in concen- 
tric, radially spaced tracks at a predetermined baud rate. When reading this recorded data, the read/write head again 
passes over the magnetic medium and transduces the magnetic transitions into pulses in an analog read signal that 
alternate in polarity These pulses are then decoded by read channel circuitry to reproduce the digital data. 

30 Decoding the pulses Into a digital sequence can be performed by a simple peak detector in a conventional analog 
read channel or. as in more recent designs, by a discrete time sequence detector in a sampled amplitude read channel. 
Discrete time sequence detectors are preferred over simple analog pulse detectors because they compensate for inter- 
symbol interference (ISI) and are less susceptible to channel noise. As a result, discrete time sequence detectors 
increase the capacity and reliability of the storage system. 

35 There are several well known discrete time sequence detection methods including discrete time pulse detection 
(DPD), partial response (PR) with Viterbi detection, maximum likelihood sequence detection (MLSD), decision-feed- 
back equalization (DFE). enhanced decision-feedback equalization (EDFE). and fixed-delay tree-search with dedsion- 
feedback (FDTS/DF). 

In conventional peak detection schemes, analog circuitry, responsive to threshold aossing or derivative Informa- 

40 tion, detects peaks in the continuous time analog signal generated by the read head. The analog read signal is "seg- 
mented" into bit cell periods and interpreted during these segments of time. The presence of a peak during the bit cell 
period is detected as a "1" bit, whereas the absence of a peak Is detected as a "0" bit. The most common errors in 
detection occur when the bit cells are not correctly aligned with the analog pulse data. Timing recovery, then, adjusts 
the bit cell periods so tiiat the peaks occur in the center of the bit cells on average in order to minimize detection errors. 

45 Since timing information is derived only when peaks are detected, the input data sti-eam is normally run length limited 
(RLL) to limit the number of consecutive "0" bits. 

As the pulses are packed closer together on the concentric data tracks in the effort to increase data density, detec- 
tion errors can also occur due to intersymbol Interference, a distortion in the read signal caused by closely spaced over- 
lapping pulses. This interference can cause a peak to shift out of Its bit cell, or its magnitude to decrease, resulting in a 

50 detection error. The ISI effect is reduced by decreasing tiie data density or by employing an encoding scheme that 
ensures a minimum number of "0" bits occur between "1" bits. For example, a (d,k) run length limited (RLL) code con- 
sti'ains to d the minimum number of "0" bits between "1" bits, and to k the maximum number of consecutive "0" bits. A 
typical (1 ,7) RLL 2/3 rate code ^codes 8 bit data words into 12 bit codewords to satisfy the (1,7) constraint. 

Sampled amplitude detection, such as partial response (PR) with Viterbi detection, allows for increased data den- 

55 sity by compensating for intersymbol interference and the effect of channel noise. Unlike conventional peak detection 
systems, sampled amplitude recording detects digital data by interpreting, at discrete time instances, tiie actual value 
of the pulse data. To this end, tiie read channel comprises a sampling 6eAce for sampling the analog read signal, and 
a timing recovery circuit for synchronizing tiie samples to the baud rate (code bit rate). Before sampling the pulses, a 
variable gain amplifier adjusts the read signal's amplitude to a nominal value, and a low pass analog filter filters the read 
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signal to attenuate aliasing noise. After sampling, a digital equalizer filter equalizes the sample values according to a 
desired partial response, and a discrete time sequence detector, such as a Viterbi detector, Interprets the equalized 
r sample values in context to determine a most likely sequence for the digital data (i.e., maximum likelihood sequence 

^ detection (MLSD)). MLSD takes into account the effect of ISI and channel noise in the detection algorithm, thereby 

5 decreasing the probability of a detection error. This increases the effective signal to noise ratio and, for a given (d,k) 
constraint, allows for significantly higher data density as compared to conventional analog peak detection read chan- 
/ nels. 

The application of sampled amplitude techniques to digital communication channels is well documented. See Y. 
Kabal and S. Pasupathy, "Partial Response Signaling", IEEE Trans. Commun. Tech., Vol. COM-23, pp.921 -934, Sept. 

10 1975; and Edward A. Lee and David G. Messerschmitt, "Digital Communication", Kluwer Academic Publishers, Boston, 
1990; and G.D. Forney, Jr., The Viterbi Algorithm". Proc. IEEE, Vol. 61 . pp. 268-278, March 1973. 

Applying sampled amplitude techniques to magnetic storage systems is also well documented. See Roy D. Cideci- 
yan. Francois Dolivo, Walter Hirt, and Wolfgang Schott. "A PRML System for Digital Magnetic Recording", IEEE Journal 
on selected Areas in Communications, Vol. 10 No. 1, January 1992, pp.38-56; and Wood et al, "Viterbi Detection of 

15 Class IV Partial Response on a Magnetic Recording Channel ", IEEE Trans. Commun., Vol. Com-34, No. 5, pp. 454- 
461. May 1986; and Coker Et al, "Implementation of PRML in a Rigid Disk Drive", IEEE Trans, on Magnetics, Vol. 27, 
No. 6, Nov. 1991;' and Carley et al. "Adaptive Continous-Time Equalization Followed By FDTS/DF Sequence Detec- 
tion", Digest of The Magnetic Recording Conference, August 15-17. 1994, pp. C3; and Moon et al, "Constrained-Com- 
plexity Equalizer Design for Fixed Delay Tree search with Decision Feedback". IEEE Trans, on Magnetics, Vol. 30, No. 

20 5, Sept. 1994; and Abbott et al, Timing Recovery For Adaptive Dedsion Feedback Equalization of The Magnetic Stor- 
age Channel", Globecom'90 IEEE Global Telecommunications Conference 1990, San Diego, CA, Nov. 1990, pp.1 794- 
1799; and Abbott et al, "Performance of Digital Magnetic Recording with Equalization and Offtrack Interference". IEEE 
Transactions on Magnetics, Vol. 27, No. 1 , Jan. 1991 ; and Cioffi et al, "Adaptive Equalization in Magnetic-Disk Storage 
Channels", IEEE Communication Magazine, Feb. 1990; and Roger Wood. "Enhanced Decision Feedback Equaliza- 

25 tion", lntermag'90. 

Similar to conventional peak detection systems, sampled amplitude detection requires timing recovery in order to 
correctly extract the digital sequence. Rather than process the continuous signal to align peaks to the center of bit cell 
periods as in peak detection systems, sampled amplitude systems synchronize the pulse samples to the baud rate. In 
conventional sampled amplitude read channels, timing recovery synchronizes a sampling clock by minimizing an error 
30 between the signal sample values and estimated sample values. A pulse detector or slicer determines the estimated 
sample values from tiie read signal samples. Even in tiie presence of ISI the sample values can be estimated and, 
together with the signal sample values, used to synchronize tiie sampling of the analog pulses in a decision-directed 
feedback system. 

A phase-Iocked-loop (PLL) normally implements the timing recovery decision-directed feedback system. The PLL 
35 comprises a phase detector for generating a phase error based on the difference between the estimated samples and 
the read signal samples. A PLL loop filter filters the phase error, and the filtered phase error operates to synchronize 
the channel samples to the baud rate. Conventionally, the phase error adjusts the frequency of a sampling clock which 
is typically the output of a variable frequency oscillator (VFO). The output of the VFO controls a sampling device, such 
as an analog-to-digital (A/D) converter, to synchronize the sampling to tiie baud rate. 
40 As mentioned above, sampled amplitude read channels also commonly employ a discrete time equalizer filter to 
equalize tiie sample values into a desired partial response (PR4, EPR4, EEPR4, etc.) before sequence detection. To 
this end. adaptive algoritiims have been applied to compensate in real time for parameter variations in the recording 
system and across the disk radius. For example, U.S. Pat. No. 5,381 ,359 entitied "Adaptation and Training of Digital 
Finite Impulse Response Filter Witiiin PRML Sampling Data Detection Channel", discloses an adaptive equalizer filter 
45 that operates accading to a well known least mean square (LMS) algoritiim. 

where W^ represents a vector of filter coefficients; ^ is a programmable gain; e^ represents a sample error between the 
50 filter's actual output and a desired output; and represents a vector of sample values from tiie filter input. In otiier 
words, the LMS adaptive equalizer filter is a closed loop feedback system that attempts to minimize the mean squared 
error between an actual output of the filter and a desired output by continuously adjusting the filter's coefficients to 
achieve an optimum frequency response. 

A problem associated with adaptive equalizer filters in sampled amplitude read channels is that the timing recovery 
55 and gain control loops can interfere with the adaptive feedback loop, thereby preventing tiie adaptive equalizer filter 
from converging to an optimal state. This non-convergence is manifested by the filter's phase and gain response drifting 
as it competes with the timing and gain control loops. An article by J.D. Coker et al. entitied "Implementation of PRML 
in a Rigid Disk Drive", published in IEEE Transactions on Magnetics, vol. 27, No. 6, Nov. 1991, suggests a three tap 
transversal filter comprising a fixed center tap and symmetric side taps in order to constrain the phase response of tiie 
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equalizer filter except in ternr^ of a fixed group delay. Constraining the phase response of the adaptive equalizer in this 
manner, however, is a very sub-optimal method for attenuating interference from the timing recovery and gain control 
loops. Furthernwre, it significantly reduces control over the adaptive filter's phase response, thereby placing the burden 
of phase compensation on the analog equalizer. 

Yet another problem associated with conventional adaptive equalizer filters is an inherent limitation on its order (i.e.. 
the number of coefficients): because the adaptive equalizer is inside the timing recovery feedback loop, its order must 
be limited to minimize the associated transport delay. Compensating for the deficiencies of the discrete time equalizer 
requires filtering the analog read signal witii a higher order analog equalizer prior to the timing recovery loop, which is 
undesirable. 

There is, tiierefore, a need for an adaptive, discrete time equalizer filter in a sampled amplitude read channel hav- 
ing an improved method for constraining the phase and gain response in order to minimize interference from the timing 
recovery and gain control loops. A further aspect of the present invention is to remove the adaptive equalizer, and its 
associated latency, from tiie timing recovery loop, tiiereby allowing a higher order discrete time filter and a simplified 
analog filter. 

SUMMARY QF TH^ INVENTION 

A sampled amplitude read channel for magnetic disk recording which asynchronously samples the analog read sig- 
nal, adaptively equalizes the resulting discrete time sample values according to a target partial response, extracts syn- 
chronous sample values through interpolated timing recovery, and detects digital data from the synchronous sample 
values using a Viterbi sequence detector is disclosed- To minimize interference from the timing and gain control loops, 
the phase and magnitude response of the adaptive equalizer filter are constrained at a predetermined frequency using 
an optimal orthogonal projection operation as a modification to a least mean square (LMS) adaptation algorithm. Fur- 
ther, with interpolated timing recovery, the equalizer filter and its associated latency are removed from the timing recov- 
ery loop, thereby allowing a higher order discrete time filter and a lower order analog filter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other aspects and advantages of the present invention will be better understood by reading the fol- 
30 lowing detailed description of the invention In conjunction witii the drawings, wherein: 

FIG. 1 is a block diagram of a conventional sampled amplitude recording channel. 

FIG. 2A shows an exemplary data format of a magnetic disk having a plurality of concentric tracks comprised of a 
plurality of user data sectors and embedded servo data sectors. 
35 FIG. 2B shows an exemplary format of a user data sector. 

FIG. 3 is a block diagram of the improved sampled amplitude read channel of the present invention comprising 
interpolated timing recovery for generating interpolated sample values and a synchronous data clock for clocking 
operation of a discrete time sequence detector and peak detector. 

FIG. 4A is a detailed block diagram of tiie prior art sampling timing recovery comprising a sampling VFO. 
40 FIG. 4B is a detailed block diagram of the interpolating timing recovery of the present invention comprising an inter- 
polator. 

FIG. 5 illustrates the channels samples in relation to the interpolated baud rate samples for the acquisition pream- 
ble. 

FIG. 6 shows an FIR filter implementation for the timing recovery interpolator. 
45 FIG. 7 depicts a cost reduced implementation for the timing recovery interpolator. 

FIG. 8A is a block diagram of a conventional adaptive, discrete time equalizer filter in a sampled amplitude read 
channel. 

FIG. 8B shows tiie adaptive, discrete time equalizer of the present invention. 

FIG. 8C shows an alternative embodiment for the adaptive, discrete time equalizer of tiie present invention. 
50 FIG. 9A illustrates the present invention adaptive filter's gain response constrained at a normalized frequency of 
1/4T 

FIG. 9B shows the present invention adaptive filter's phase response constrained at a normalized frequency of 
1/4T 

FIG. 1 0 illustrates operation of an orthogonal projection operation of the present invention for constraining the gain 
55 and phase response of tiie adaptive filter 

FIG. 1 1 A shows an implementation for a reduced cost orthogonal projection operation. 
FIG. 1 1 B shows an alternative embodiment for the reduced cost orthogonal projection operation. 
FIG. 1 1C illustrates an implementation for a gradient averaging circuit used in the reduced cost orthogonal projec- 
tion operation of FIG. 1 1 B. 
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FIG. 1 1D shows yet another alternative embodiment for the reduced cost orthogonal projection operation of frie 
present invention. 

DgTAIL^D DESCRIPTION QF THE PREFERRED EMBODIMENT 

5 

Conventional Sampled Amplitude Read Channel 

Referring now to FIG. 1 , shown is a detailed block diagram of a conventional sanpled amplitude read channel. Dur- 
ing a write operation, either user data 2 or preamble data from a data generator 4 (for example 2T preamble data) is 

10 written onto the media. An RLL encoder 6 encodes the user data 2 into a binary sequence b(n) 8 according to an RLL 
constraint. A precoder 10 precodes the binary sequence b{n) 8 in order to compensate for the transfer function of the 
recording channel 18 and equalizer filters to form a preceded sequence -b(n) 12. The preceded sequence -b(n) 12 is 
converted into symbols a(n) 16 by translating 14 ~b(N) = 0 into a(N) = -1, and -b^N) = 1 into a(N) = +1. Write circuitry 
9, responsive to the symbols a(n) 1 6, modulates the cun^ent in the recording head coil at the baud rate 1/T to record the 

15 binary sequence onto the media. A frequency synthesizer 52 provides a baud rate write clock 54 to the write circuitry 9 
and is adjusted by a channel data rate signal (CDR) 30 according to the zone the recording head is over. 

When reading the recorded binary sequence from the media, timing recovery 28 first locks to the write frequency 
by selecting, as the input to the read channel, the write clock 54 through a multiplexor 60. Once locked to the write fre- 
quency, the multiplexor 60 selects the signal 19 from the read head as the Input to the read channel in order to acquire 

20 an acquisition preamble recorded on the disk prior to the recorded user data. A variable gain amplifier 22.adjusts the 
amplitude of the analog read signal 58, and an analog filter 20 provides initial equalization toward the desired response 
as well as attenuating aliasing noise. A sampling device 24 samples the analog read signal 62 from the analog filter 20. 
and a discrete time equalizer filter 26 provides further equalization of the sample values 25 toward the desired 
response. In partial response recording, for example, the desired response is often selected from Table 1 . 

25 The discrete equalizer filter 26 may be implemented as a real-time adaptive filter which compensates for parameter 
variations over the disk radius (i.e., zones), disk angle, and environmental conditions such as temperature drift. To this 
end, the filter 26 receives estimated sample values B143 generated by tiie timing recovery circuit 28; the estimated 
samples being input into an adaptive feedback loop and used to generate sample en-ors. The adaptive feedback loop 
conventionally employs a least mean square (LMS) algorithm to adapt the filter coefficients (i.e.. it adapts the frequency 

30 and phase response of the filter) until a minimum sample error is achieved. Operation of a conventional adaptive equal- 
izer filter is discussed in greater detail below. 

After equalization, the equalized sample values 32 are applied to a decision directed gain control 50 and timing 
recovery 28 circuit for adjusting tine amplitude of the read signal 58 and the frequency and phase of the sampling device 
24, respectively. Timing recovery adjusts the frequency of sampling device 24 over line 23 in order to synchronize tiie 

35 equalized samples 32 to the baud rate. Frequency synthesizer 52 provides a course center frequency setting to the tim- 
ing recovery circuit 28 over line 64 in order to center tiie timing recovery frequency over temperature, voltage, and proc- 
ess variations. The channel data rate (CDR) 30 signal adjusts a frequency range of the synthesizer 52 according to the 
data rate for the current zone. Gain control 50 adjusts the gain of variable gain amplifier 22 over line 21 in order to match 
the magnitude of the channel's frequency response to the desired partial response. 

40 The equalized samples Y(n) 32 are also sent to a discrete time sequence detector 34, such as a maximum likeli- 
hood (ML) Viterbi sequence detector, which detects an estimated binary sequence *b(n) 33 from the sample values. An 
RLL decoder 36 decodes the estimated binary sequence '^(n) 33 from the sequence detector 34 into estimated user 
data 37. A data sync detector 66 detects the sync mark 70 (shown in FIG. 2B) in the data sector 15 in order to frame 
operation of the RLL decoder 36. In tiie absence of errors, the estimated binary sequence *b(n) 33 matches tiie 

45 recorded binary sequence b(n) 8. and the decoded user data 37 matches tiie recorded user data 2. 

Data Format 

FIG. 2A shows an exemplary data format of a magnetic media comprising a series of concentric data tracks 13 
50 wherein each data track 13 comprises a plurality of sectors 15 with embedded servo wedges 17. A servo confroller (not 
shown) processes the servo data in the servo wedges 17 and, in response thereto, positions the read/write head over 
a desired track. Additionally, tiie servo controller processes servo bursts within tiie servo wedges 17 to keep the head 
aligned over a centerline of tiie desired track while writing and reading data. The servo wedges 1 7 may be detected by 
a simple discrete time pulse detector or by the discrete time sequence detector 34. If tiie sequence detector 34 detects 
55 the servo data, then the format of the servo wedges 1 7 includes a preamble and a sync marK similar to the user data 
sectors 15. 

FIG. 2B shows tiie format of a user data sector 15 comprising an acquisition preamble 68, a sync mark 70, and 
user data 72. Timing recovery uses tiie acquisition preamble 68 to acquire the correct sampling frequency and phase 
before reading tiie user data 72, and the sync mark 70 demarks the beginning of the user data 72 (see co-pending U.S. 
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patent application 08/313,491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording!. 

To increase the overall storage density, the disk Is partitioned Into an outer zone 1 1 comprising fourteen data sec- 
tors per track, and an inner zone 27 comprising seven data sectors per track. In practice, the disk is actually partitioned 
Into several zones with a different number of sectors In each zone, and the data recorded and detected at a different 
5 data rate in each zone. 

Improved Sampled Amplitude Read Channel 

FIG. 3 shows the inproved sampled amplitude read channel of the present invention wherein the conventional 
10 sampled timing recovery 28 of FIG. 1 has been replaced by interpolated timing recovery B100. In addition, the write 
frequency synthesizer 52 generates a baud rate write clock 54 applied to the write circuitry 9. and an asynchronous 
read clock 54 for clocking the sampling device 24, the discrete time equalizer filter 26. artd the interpolated timing recov- 
ery B100 at a frequency relative to the current zone (CDR 30). In an alternative embodiment, a first frequency synthe- 
sizer generates the write clocK and a second frequency synthesizer generates the read clock. 
IS The discrete equalizer filter B103 is real-time adaptive, receiving Interpolated sample values 81 02 and estimated 
sample values 8143 from tiie interpolated timing recovery circuit 8100 for use in a modified least mean square (LMS) 
algorithm which constrains the filter's gain and phase response according to the present invention, the details of which 
are set forth below. 

The interpolated timing recovery 8100 interpolates tiie equalized sample values 32 to generate interpolated sam- 
20 pie values 8102 substantially synchronized to the data rate of the current zone. A discrete time sequence detector 34 
detects an estimated binary sequence 33 representing the user data from the interpolated sample values 8102. The 
interpolated timing recovery 8100 circuit generates a synchronous data dock 8104 for clocking operation of the gain 
control 50. discrete time sequence detector 34. sync mark detector 66 and RLL decoder 36. 

25 Conventional Timing Recovery 

An overview of tiie conventional sampling timing recovery 28 of FIG. 1 is shown in FIG. 4A. The output 23 of a var- 
iable frequency oscillator (VFO) 81 64 conti-ols the sampling clock of a sampling device 24 which is typically an analog- 
to-digital converter (A/D) in digital read channels. A multiplexor 81 59 selects tiie unequalized sample values 25 during 

30 acquisition and tiie equalized sample values 32 during tracking, thereby removing the discrete equalizer filter 26 from 
the timing loop during acquisition in order to avoid its associated latency. A phase error detector 8155 generates a 
phase error in response to tiie sample values received over line 8149 and estimated sample values -^Y^ from a sample 
value estimator 8141. such as a slicer in a d=0 PR4 read channel, over line 8143. A loop filter 8160 filters the phase 
error to generate a frequency offset Af 8167 tiiat setties to a value proportional to a frequency difference between the 

35 sampling clock 23 and the baud rate. The frequency offset Af 8167, together with the center frequency control signal 
64 from the frequency synthesizer 52. adjust the sampling clock 23 at the output of the VFO 8164 in order to synchro- 
nize the sampling to the baud rate. 

A zero phase start 81 62 circuit suspends operation of the VFO 81 64 at the beginning of acquisition in order to min- 
imize the initial phase error between the sampling clock 23 and the read signal 62. This is achieved by disabling the 

40 VFO 81 64, detecting a zero crossing in the analog read signal 62, and re-enabling the VF0 1 64 after a predetermined 
delay between ttie detected zero crossing and tine first baud rate sample. 

The estimated sample values 8143 at The output of the slicer 8141 are also input into the discrete time equalizer 
filter 26 of FIG. 1 for use in a conventional least mean square (UVIS) adaptation algorithm as is described in more detail 
below. 

45 

Interpolated Timing Recovery 

The interpolated timing recovery 8100 of tiie present invention is shown in FIG. 48. The VFO 8164 in the conven- 
tional timing recovery of FIG. 4A is replaced with a modulo-Ts accumulator 8120 and an interpolator 8122. In addition, 

50 an expected sample value generator 81 51 , responsive to interpolated sample values 81 02. generates expected sam- 
ples used by the phase error detector 8155 to compute the phase enror during acquisition. A multiplexor 8153 
selects tiie estimated sample values -Yk+^ from the slicer 8141 for use by the phase error detector 8155 during track- 
ing. The data clock 81 04 is generated at the output of an AND gate 8126 in response to the sampling clock 54 and a 
mask signal 81 24 from the modulo-Ts accumulator 81 20 as discussed in further detail below. The phase error detector 

55 8155 and tiie slicer 8141 process interpolated sample values 8102 at the output of the interpolator 8122 rather tiian 
the channel sample values 32 at the output of the disaete equalizer filter 26 as in FIG. 4A. A PID loop filter 8161 con- 
trols the closed loop frequency response similar to the loop filter 8160 of FIG. 4A. 

The interpolated sample values Yk+t B102 and the estimated sample values -Yk+T from the slicer 8141 are input 
into the adaptive, discrete equalizer filter 8103 of FIG. 3 for use by a modified least mean square (l-MS) algorithm, the 
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details of which are set forth below. 

In the interpolated timing recovery of the present invention, locking a VFO to a reference frequency before acquir- 
ing the preanrtble is no longer necessary; multiplexing 60 the write dock 54 into the analog receive filter 20 (as In FIG. 
1) is not necessary. Further, the sampling device 24 and the discrete equalizer filter 26. together with their associated 

5 delays, have been removed from the timing recovery loop; it is not necessary to multiplex B159 around the equalizer 
filter 26 between acquisition and tracking. However, it is still necessary to acquire a preanrtble 68 before tracking the 
user data 72. To this end, a zero phase start circuit B163 minimizes the initial phase error between the interpolated 
sanple values and the baud rate at the beginning of acquisition similar to the zero phase start circuit B162 of FIG. 4A. 
However, rather than suspend operation of a sampling VFO 81 64. the zero phase start circuit 81 63 for interpolated tim- 

10 ing recovery computes an initial phase error x from the A/D 24 sample values 25 and loads this initial phase error into 
the modulo-Ts accumulator 8120. 

For more details concerning the P ID loop filter 81 61 , phase error detector 81 55. expected sample generator 81 51 , 
and slicer 8141 , refer to the above referenced co-pending U.S. patent applications "Sampled Amplitude Read Channel 
Comprising Sample Estimation Equalization, Defect Scanning, Channel Quality. Digital Servo Demodulation, PID Filter 

15 for Timing Recovery, and DC Offset Control" and "Improved Timing Recovery For Synchronous Partial Respor^e 
Recording." A detailed description of the modulo-Ts accumulator 81 20, data clock 81 04, and interpolator 81 22 is pro- 
vided in tiie following discussion. 

Interpolator 

20 

The interpolator 81 22 of FIG. 4B is understood with reference to FIG. 5 which shows a sampled 2T acquisition pre- 
amble signal 8200. The target synchronous sanple values 8102 are shown as black circles and the asynchronous 
channel sample values 32 as vertical arrows. Beneath the sampled preamble signal is a timing diagram depicting tiie 
corresponding timing signals for the sampling clock 54, the data clock 81 04 and the mask signal 81 24. As can be seen 
25 in FIG. 5. the preamble signal 8200 is sampled slightiy faster than tiie baud rate (the rate of the target values). 

The function of the interpolator is to estimate the target sample value by interpolating the channel sample values. 
For illustrative purposes, consider a simple estimation algorithm, linear interpolation: 

Y(N-1) = x(N-1) + T • (x(N) - x(N-1)); where: (1) 

30 

x(N-1) and x(N) are tiie channel samples sun'ounding the target sample; and x is an interpolation interval proportional 
to a time difference between the channel sample value x(N-1) and the target sample value. The interpolation interval x 
is generated at the output of modulo-Ts accumulator 8120 which accumulates the frequency offset signal Af 8167 at 
the output of the PID loop filter 8161 : 

35 

X = (ZAO MOD TS; where: (2) 

Ts is the sampling period of tiie sanrrpling clock 54. Since tiie sampling clock 54 samples the analog read signal 62 
slightly faster than the baud rate, it is necessary to mask the data clock every time the accumulated frequency offset Af, 

40 integer divided by Ts, increments by 1. Operation of tiie data clock 8104 and tiie mask signal 8124 generated by tiie 
modulo-Ts accumulator 8120 is understood with reference to the timing diagram of FIG. 5. 

Assuming the interpolator implements the simple linear equation (1) above, then channel sample values 8202 and 
8204 are used to generate the interpolated sample value corresponding to target sample value 8206. The interpolation 
interval x 8208 is generated according to equation (2) above. The next interpolated sample value corresponding to tiie 

45 next target value 821 0 is computed from channel sample values 8204 and 821 2. This process continues until tiie inter- 
polation interval x 8214 would be greater than Ts except that it "wraps" around and is actually x 821 6 (i.e., the accumu- 
lated frequency offset Af, integer divided by Ts, increments by 1 causing the mask signal 81 24 to activate). At tiiis point, 
the data clock 8104 is masked by mask signal 8124 so tiiat tiie interpolated sample value con^esponding to the target 
sample value 8220 is computed from channel sample values 8222 and 8224 ratiier tiian channel sample values 8218 

50 and 8222. 

The simple linear interpolation of equation (1) will only work if the analog read signal is sampled at a much higher 
frequency than tiie baud rate. This is not desirable since operating tiie channel at higher frequencies increases its com- 
plexity and cost. Therefore, in tiie preferred embodiment the interpolator 8122 is implemented as a filter responsive to 
more than two channel samples to compute the interpolated sample value. 
55 The ideal discrete time phase interpolation filter has a flat magnitude response and a constant group delay of x: 

C,(eh = e^"' (3) 

which has an ideal impulse response: 
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sinc(7i-(n-xyTg)). (4) 

Unfortunately, the above non-causal infinite impulse response (4) cannot be realized. Therefore, the impulse response 
of the interpolation filter is designed to be a best fit approximation of the ideal Impulse response (4). This can be accom- 
5 plished by minimizing a mean squared en-or between the frequency response of the actual interpolation filter and the 
frequency response of the ideal interpolation filter (3). This approximation can be improved by taking into account the 
spectrum of the input signal, that is. by minimizing the mean squared error between the input spectrum multiplied by 
the actual interpolation spectrum and the Input spectrum multiplied by the ideal Interpolation spectrum: 

10 C,(ehx(eh - C,(e^'')X(eh: where: (5) 

C^(ej'°) Is the spectrum of the actual interpolation filter; and X(e^ is the spectrum of the Input signal. From equation 
(5), the mean squared error is represented by: 

15 n 

-n 



20 X{e^^) Is the spectrum of tiie read channel (e.g., PR4, EPR4, EEPR4 of Table 1 or some other partial response spec- 
trum). 

In practice, the above mean squared error equation (6) is modified by specifying tiiat the spectrum of tiie Input sig- 
nal is bandlimitted to some predetermined constant O^oa^an where 0 < a < 1 ; that is: 

25 |X(ehl = 0.for|Q)|^an. 

Then equation (6) can be expressed as: 

an 
-an 



The solution to the minimization problem of equation (7) Involves expressing the actual interpolation filter in terms of its 
35 coefficients and then solving for tiie coefficients that minimize the error in a classical mean-square sense. 
The actual Interpolation filter can be expressed as the FIR polynomial: 

n=R-1 

5^(6^"')= X '^"^^ : where: (8) 

40 n=-R 



2 R is the number of taps in each interpolation filter and the sample period Ts has been normalized to 1 . A mathematical 
derivation for an interpolation filter having an even number of coefficients is provided below. It is wltiiin the ability of 
45 those skilled In the art to modify the mathematics to derive an interpolation filter having an odd number of coefficients. 
Substituting equation (8) into equation (7) leads to the desired expression In terms of tiie coefficients C^(n): 

an rNR-1 

E'a = 2b 1 I S Gx(n)e"^^" - e^-'^I^IXCehl^c/co. (9) 

-an n=-R 

The next step is to take the derivatives of equation (9) with respect to the coefficients C^(n) and set them to zero: 

-Oforno = -R. ..,,0. 1..... R-1. (10) 



ac,(no) 

After careful manipulation, equation (10) leads to: 
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ait pxn«R-1 N -| 

J Z C»cos(a)(n^ - n)) - cos( <o{n^ + t)) IXCe^l'do) = 0 



(11) 



for no = -R 0, 1 R-1 

Defining (|>(r) as: 



10 



m= J |X(enr cos(a)r) tf© 



(12) 



15 



20 



25 



30 



35 



40 



and substituting equation (12) into equation (1 1) gives: 



n«R-1 

X C^(n)(|)(n-nJ = (l,(n^ + T) 

n=-R 



(13) 



for no = -R. .... 0. 1, .... R-1 . 

Equation (13) defines a set of 2R linear equations in terms of the coefficients C^(n). Equation (13) can be expressed 
more compactly in matrix form: 



is a column vector of the form: 



<1>T is a Toeplitz matrix of the form: 



<t>jC^ = 4),; where: 



C, = [c,(-R) c,(0) c,(R-1)]' 



(t)(0) (^(1) 
0(1) (|>(0) 

<j»(2R-l) 



... <1)(2R-1 



.. (1)(0) 



and is a column vector of the form: 
45 = [(K-R + x) <t>(T). (|>(1 + x) (t»(R-1 + x)l \ 



50 



The solution to equation (14) is: 



= ^j'^O^; where: 



(14) 



(15) 



0f ^ is an inverse matrix that can be solved using well known methods. 



Table B2 shows example coefficients C^{n) calculated from equation (15) with 2R=6. a=0.8 and X(e^ )=PR4 . The 
implementation of the six tap FIR fitter is shown In FIG. 6. A shift register B250 receives the channel samples 32 at the 
sanpllng clock rate 54. The filter coefficients C^(n) are stored in a coefficient register file B252 and applied to cone- 
55 sponding multipliers according to the current value of x B128. The coefficients are multiplied by the channel samples 
32 stored in the shift register 8250. The resulting products are summed 8254 and the sum stored in a delay register 
8256. The coefficient register file 8252 and the delay register 8256 are clocked by the data clock 8104 to implement 
the masking function described above. 

In an alternative embodiment not shown, a plurality of static FIR filters, having coefficients that correspond to the 
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different values of x, filter the sanple values In the shift register B250. Each filter outputs an interpolation value, and the 
cun-ent value of the Interpolation interval x B1 28 selects the output of the corresponding f ilt©' as the output B1 02 of the 
interpolator 81 22. Since the coefficients of one filter are not constantly updated as in FIG. 6, this multiple filter embod- 
iment increases the speed of the Interpolator 8122 and the overall throughput of the read channel. 

5 

Cost Reduced Interpolator 

Rather than store all of the coefficients of the interpolation filters in memory, in a more efficient, cost reduced imple- 
mentation the coefficient register file 8252 of FIG. 6 computes the filter coefficients C^(n) in real time as a function of x. 
10 For example, the filter coefficients C^(n) can be computed in real time according to a predetermined polynomial in x 
(see, for example, U.S. Pat. No. 4,866.647 issued to Farrow entitled. "A Continuously Variable Digital Delay Circuit," the 
disclosure of which is hereby incorporated by reference). An alternative, preferred embodiment for computing the filter 
coefficients in real time estimates the filter coefficients according to a reduced rank matrix representation of the coeffi- 
cients. 

15 The b»ank of filter coefficients stored in the coefficient register file 8252 can be represented as an MxN matrix A^^xN* 
where N is the depth of the interpolation filter (i.e.. the number of coefficients C^(n) in the impulse response computed 
according to equation (15)) and M is the number of interpolation Intervals (i.e., the number of x intervals). Rather than 
store tiie entire Ai^^n matrix in memory, a more efficient, cost reduced implementation is attained through factorization 
and singular value decomposition (SVD) of the A^xn matrix. 

20 Consider that the A^xn matrix can be factored into an Fpyix^ and Gnxn matrix. 

^MxN = ^MxN * ^NxN • 

Then a reduced rank approximation of the A^xn matrix can be formed by reducing the size of the F^xn and Gnxn matri- 
25 ces by replacing N with L where L<N and, preferably, L«N. Stated differently, find the F^xl and G|_xn matrices whose 
product best approximates the A|^xn "T^iatrix, 

^MxN * ^MxL * ^LxN • 

30 The convolution process of the Interpolation filter can then be carried out. as shown in FIG. 7, by implementing tiie Glxn 
matrix as a bank of FIR filters 8260 connected to receive the channel sample values 32. and the F^^l matrix imple- 
mented as a lookup table 8262 indexed by x 8128 (as will become more apparent in the following discussion). Those 
skilled in the art will recognize that. In an alternative embodiment, the Ai^^n matrix can be factored Into more than two 
matrices (i.e.. A « FGH...). 

35 The preferred method for finding tiie F^xL and G|_xn matrices is to minimize the following sum of squared en'ors: 

M N 

I I(Ajn-(FMxL •GLxN)]n)' 06) 

j= 1n=1 

40 

The solution to equation (16) can be derived through a singular value decomposition of the A^xn matrix, comprising tiie 
steps of: 

45 1 . performing an SVD on the A^xn matrix which gives tiie following unique factorization (assuming M^N): 

A MxN = U MxN • D NxN • V NxN where: 

Umxn is a f^xN unitary matrix; 

50 Dnxn 'S a NxN diagonal matrix (cti , aN) where aj are tiie singular values of A|^xn» and 02-^0^'^ 

0; and 

Vnxn 'S a NxN unitary matrix; 

2. selecting a predetermined L number of tiie largest singular values a to generate a reduced size diagonal matrix 
55 D|_xl: 
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10 



Dlxl =Diag {apGj, ...,aL} = 



0 .. 
0 ^2 0. 



0 



0 



. 0 

0 Or 



15 



20 



25 



3. extracting the first L columns from the U^xN matrix to form a reduced U^xl matrix: 



U 



MxL - 



... u 



1,N 



U 



M,N 



30 



35 



4. extracting the first L rows from the V^xn "matrix to form a reduced Vl^n matrix: 



^LxN - 



40 




45 

5. defining the F^xl and Glxn matrices such that: 

^ MxL * ^ LxN = MxL * ^ LxL ' ^ UN ^ MxN 

50 (for example, let F = U * ^ ^l and G lxn = V lxn )- 

In the above cost reduced polynomial and reduced rank matrix embodiments, the interpolation filter coefficients 
C^(n) are computed in real time as a function of t; that is. the filter's impulse response h(n) is approximated according to: 

55 l 

h(n,t) = c^(n)= ^G,(n) • f(l,x); where: (17) 
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f(i,T) is a predetermined furKition in x (e.g., polynomial in x or x indexes the above F^xl nnatrlx); L is a degree which 
determines the accuracy of the approximation (e.g., the order of the polynomial or the column size of the above F^xt 
matrix); and Gi(n) is a predetermined matrix (e.g.. the coefficients of the polynomial or the above Gijdsi matrix). As L 
increases, the approximated filter coefficients C^(n) of equation (17) tend toward the ideal coefficients derived from 
5 equation (15). It follows from equation (17) that the output of the interpolation f Oter Y(x) can be represented as: 

N L 

Y(x)= U(x-n) 5;G,(n) • f(i.x) (18) 
ns1 i=1 

10 

where U(x) are the channel sample values 32 and N is the number of interpolation filter coefficients C^(n). 

Referring again to FIG. 6. the coefficient register file can compute the interpolation filter coefficients C^(n) according 
to equation (17) and then convolve the coefficients C^(n) with the channel samples U(x) 32 to generate the interpolated 
sample values B102 synchronized to the baud rate. However, a more efficient Implementation of the interpolation fitter 
15 can be achieved by rearranging equation (1 8) : 

L N 

Y(x)= XKi.^) ' 

1=1 n= 1 

20 

FIG. 7 shows the preferred embodiment of the interpolation filter according to equation (19). In the polynomial 
embodiment, the function of x is a polynomial in x, and the matrix Gj(n) are the coefficients of the polynomial. And in the 
reduced rank matrix embodiment, the function of x is to index the above F^^l n^atrix B262, and the second summation 
25 in equation (19), 

N 

^Gj(n) • U(x-n) 

n = 1 

30 

is implemented as a bank of FIR filters B260 as shown in FIG. 7. Again, in equation (19) L is the depth of the approxi- 
mation function f(i,x) (e.g., the order of the polynomial or the column size of the above F^^l "matrix) and N is the depth 
of the interpolation filter's impulse response (i.e., the number of coefficients in the impulse response). It has been deter- 
35 mined that N=8 and L=3 provides the best performance/cost balance; however, these values may increase as IC tech- 
nology progresses and the cost per gate decreases. 

Conventional Adaptive Equalizer 

40 FIG. 8 A illustrates a prior art adaptive, discrete time equalizer that operates according to the well known least mean 
square (LMS) algorithm, 

WK.i=W,.n.eK-X,. 

45 or alternatively. 

Wk,i=WK-fi-X,.eK 

where represents a vector of FIR filter coefficients; ^ is a programmable gain; e^ represents a sample error (or vec- 
50 tor of sample errors ejj between the FIR filter's actual output and a desired output; and represents a vector of sam- 
ples values (or a scalar X^J from the FIR fitter input. To better understand operation of the present invention, the second 
representation of the I^S algorithm is used throughout this disclosure. 

TTie desired filter output is the estimated sample values -Yj^ at the output of slicer B141 . The estimated sample val- 
ues ~Yk are subtracted from the FIR fitter's output Y^ to generate the sample error ek- The UVIS algorithm attempts to 
55 minimize the sample error in a least mean square sense by adapting the FIR filter coefficients; that is, it adjusts the FIR 
filter's gain and phase response so that the overall channel response adequately tracks the desired partial response 
(e.g., PR4. EPR4, EEPR4. etc.). 

As previously mentioned, interference from the timing recovery 28 and gain control 50 loops can prevent the adap- 
tive, discrete time equalizer 26 from converging to an optimal state. For example, a phase adjustment in the adaptive 
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filter 26 can affect the sampling phase enror for timing recovery 28. Timing recovery 28 compensates for the filter's 
phase adjustment by adjusting its sampling phase: this adjustment can result in yet another phase adjustment by the 
adaptive equalizer 26. Thus, the phase response of the adaptive equalizer may never converge. Similarly, the gain con- 
trol loop 50 can interfere with the gain response of the adaptive filter 26 and prevent it from converging. 

5 

Constrained Adaptive Equalizer 

FIG. 8B illustrates operation of the adaptive, discrete time equalizer filter B103 of the present invention. Sample val- 
ues from the A/D 24 are input over line 25 to a discrete time FIR fitter CI 00 comprising a predetermined number of coef- 

10 ficients. the values of which determine the filter's gain and phase response. Because the FIR filter CI 00 operates on 
the sample values prior to the interpolated timing recovery loop B1 00, its order can be increased over the prior art with- 
out adversely affecting the latency of timing recovery (i.e., the number of filter coefficients can be Increased). 

The output 32 of the FIR filter CI 00 is input into the Interpolator B122 for generating the interpolated sample 
values 8102. The interpolated sample values Y^^.^ 81 02 are input into a slicer 8141 (FIG. 4B) which generates 

15 estimated sample values -Y^^^. The estimated sample values -Y^^^ are subtracted from tfie interpolated sample val- 
ues Y)^^^ at adder CI 02 to generate a sample error value en+x CI 04 that is synchronized to the baud rate rather than 
the sample rate. Because the LMS algorithm operates on sample values at the sample rate, it is necessary to con- 
vert tiie en-or value ek+^ CI 04 into an error value ek C112 synchronous to the sample rate. This is accomplished by an 
interpolation circuit CI 06 which computes an interpolated error value ek C112 from the baud rate error values Ck+x 

20 CI 04. Preferably, the en'or value interpolation circuit CI 06 is implemented as a first order linear interpolation, but it may 
be a simple zero order hold, or a more complex interpolation filter as described above. 

FIG. 80 shows an alternative embodiment for generating the error value ek- As illustrated, tiie estimated ideal sam- 
ple values '-Yk+^ from the slicer B141 are interpolated by the interpolation circuit CI 06 to generate estimated ideal sam- 
ple values -Yk CI 20 which are subti-acted CI 02 from the equalized sample values Yk 32 at the output of tiie FIR filter 

25 CI 00 to generate the error value ek 

In both embodiments, the error value ek C1 12 is input into a modified LMS circuit C1 14 which computes updated 
filter coefficients Wk C116 according to, 

WK.i=WK-n-Pviv/-(X,-eK) 

30 

where the operation Pv^ V2-'- is an orthogonal projection operation which constrains the gain and phase response of the 
FIR filter CI 00 in order to attenuate interference from the gain and timing loops. 

Operation of the orthogonal projection operation Py^\/2^ will now be described in relation to the gain and phase 
response of the FIR filter CI 00. FIG.s 9 A and 9B show an example gain and phase response, respectively, for the FIR 
35 filter CI 00 of the present invention. The gain and phase response vary over time as the filter adapts to parameter 
changes in the recording system; that is, the filter continuously adjusts the channel's overall frequency response so that 
it matches the desired partial response (PR4. EPR4, EEPR4, etc.) as best possible. In the present invention, interfer- 
ence from the timing recovery and gain control loops is attenuated by consti-aining the gain and phase response of tiie 
FIR filter CI 00 at a predetermined frequency using an ortiiogonal projection operation Pv^vg-*-. 
40 Referring again to FIG. 9A, the gain (magnitude) response of the FIR filter CI 00 has been constrained to a prede- 
termined value (denoted by g) at the normalized frequency of 0.5 (1/4Ts). Similarly, tiie phase response of the FIRfilter 
CI 00 has been consti-ained to kn at the normalized frequency of 0.5 as shown in FIG. 9B. In effect, these constraints 
allow tiie gain and phase to vary (adapt) at all frequencies except at the normalized frequency of 0.5. thereby constrain- 
ing the filter's frequency response in a manner that attenuates interference from the gain and timing loops. The gain 
45 constraint g is relatively arbitrary except that it is selected to optimize tiie dynamic range of tiie filter's coefficients. How- 
ever, constraining tiie filter's response at tiie normalized frequency of 0.5 and selecting a phase constraint of kir 
reduces tiie complexity of tiie orthogonal projection operation Pv^ vg-*-. and simplifies implementation of the zero phase 
start circuit B163 (FIG. 4B) and sync mark detector 66 (FIG. 1). 

As mentioned above, the zero phase start circuit B163 (FIG. 4B) minimizes tiie initial phase error between the inter- 
fa polated sample values and the baud rate at the beginning of acquisition by computing an initial phase en'or x from the 
A/D 24 sample values 25 and then loading tiiis initial phase error into tiie modulo-Ts accumulator B120. To compute 
the initial phase error x, tiie zero phase start circuit 8163 must take into account tiie phase delay of tiie adaptive equal- 
izer filter 81 03 since interpolated timing recovery 81 00 operates on the equalized samples 32. not the A/D 24 samples 
25. With the acquisition preamble 68 (FIG. 2B) having a frequency of 1/4T (i.e., 0.5 normalized), constraining tiie phase 
55 response of tiie adaptive equalizer 8103 at tiie preamble frequency (1/4Ts) fixes the phase delay of tiie equalizer B103 
during acquisition, thereby allowing tiie zero phase start circuit 8163 to accurately compute the initial phase error x. 
Furthermore, since tiie phase constraint at the preamble frequency is fixed at kjc, tiie phase delay of tiie equalizer 81 03 
will either be zero or 180^ (i.e., the adjustment to the initial phase error is nothing or a sign change). 

Constraining the phase response of tiie adaptive equalizer 8103 to kn at the preamble frequency also simplifies 



13 



EP0 805 447 A2 

irrplementation of the sync mark detector 66 (FIG. 1) in sampled amplitude read channels that use the conventional 
synchronous sampling timing recovery 28 (FIG. 4A) rather than interpolated timing recovery B100 (FIG. 4B). Operation 
of the sync mark detector 66 is described in detail in the above referenced co-pending U.S. Patent Application 
08^3.797 entitled "Improved Fault Tolerant Sync Mark Detector For sampled Amplitude Magnetic Recording." As 

5 desaibed therein, the sync mark detector 66 is enabled coherent with the end of the acquisition preamble 68 and rela- 
tive to the transport delay from the output of the A/D 24 to the sync mark detector 66. With an adaptive equalizer, the 
transport delay will vary unless the filter's phase response is constrained at the acquisition preamble frequency by, for 
example, using an orthogonal projection operation Py^y2^ of the present invention. 

Turning now to the implementation details of the orthogonal projection operation PviV2-^, the equalizer's frequency 

10 response is 



C(e''"b = XCke''^"'' 

k 

15 

where are the coefficients of the equalizer's impulse response. At the preamble frequency (1/4T), the equalizer's fre- 
quency response is 



25 where the sampling period has been normalized to T=1 . In matrix form, the equalizer's frequency response at the pre- 
amble frequency is, 



30 



35 



e 



40 



( -"AO 
e 

V J 



.nv 
^2 



(N-l) 



(j)° " 

U)-' 




1 

-j 
-1 


.0) 




j 



Those skilled in the art will recognize that shifting the time base will lead to four different, but functionally equivalent, 
frequency responses at the preamble frequency (i.e., [1,-j.-1,j,...]C, [-j,-1,j,1,...l£ [-1,j,1,-j,-.]C and D.IH.-''. ■)£)■ Con- 
45 straining the phase response of the equalizer B103 to an integer multiple of % at tiie preamble frequency (1/4T) implies 
that the imaginary component of its frequency response is zero, 
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55 



0 

-1 

0 

1 



= C^Vj =0 



If the imaginary component of the frequency response is constrained to zero, as described above, then constraining the 
magnitude of the equalizer to g at the preamble frequency (1/4T) implies tiiat the real component of the frequency 
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response equals g, 



1 
0 

-1 
0 



Therefore, the equalizer's coefficients must be constrained to satisfy the following two conditions: 

C'^-Vi =0andC'^-V2 = g. 

The above constraints are achieved by multiplying the computed gradient Xk • ^ by an orthogonal projection operation 
PviV2-^ as part of a modified LMS algorithm C114. 

To understand the operation of the orthogonal projection operation, consider an equalizer that comprises only two 
coefficients: Cq and as shown in FIG. 10. The phase consti*aint condition • V ^ =0 implies that the filter coeffi- 
cient vector must be orthogonal to . When using an unmodified LMS algorithm to update tiie filter coefficients, the 
orthogonal constraint is not always satisfied as shown in FIG. 10. The present invention, however, consti'ains the fitter 
coefficients to a subspace <£) which is ortiiogonal to by multiplying the gradient values • by a projection oper- 
ation Pvi"'-, where the null space of the projection operation Pv/ is orthogonal to <C>. The updated coefficients core- 
spond to a point on tiie orthogonal subspace <C ) closest to the coefficients derived from the unmodified LMS algorithm 
as shown in FIG. 10. 

Similar to tiie phase constraint projection operation Pv-|-*-, a second ortiiogonal projection operation Pv2-^ constrains 
the filter coefficients such tiiat the coefficient vector C"*^ satisfies the above gain const-aint: C ^ • V 2 = g . The combined 
orthogonal projection operation Pv-^V2^ eliminates two degrees of freedom in an N-dimensional subspace where N is 
the number of filter coefficients (i.e., tiie orthogonal projection operation PviV2-^ has a rank of N-2). 

An orthogonal projection operation for and V2 can be computed according to 



Pv/ = I-Pv, = l-V,(v/V,)-V/ 



(20) 



where Pv^Vg ''"=Pvi "'" • Pvg since V-| is orthogonal to V2. The orthogonal projection operation PviV2-^ computed 
35 using the above equation for an equalizer comprising ten filter coefficients is a matrix 



40 



45 



50 



55 



4 


0 


1 


0 


-1 


0 


1 


0 


-1 


0 


0 


4 


0 


1 


0 


-1 


0 


1 


0 


-1 


1 


0 


4 


0 


1 


0 


-1 


0 


1 


0 


0 


1 


0 


4 


0 


1 


0 


-1 


0 


1 


= -1 


0 


1 


0 


4 


0 


1 


0 


-1 


0 


0 


-1 


0 


1 


0 


4 


0 


1 


0 


-1 


1 


0 


-1 


0 


1 


0 


4 


0 


1 


0 


0 


1 


0 


-1 


0 


1 


0 


4 


0 


1 


-1 


0 


1 


0 


-1 


0 


1 


0 


4 


0 


0 


-1 


0 


1 


0 


-1 


0 


1 


0 


4 



The above matrix Pv^V2-'- is an orthogonal projection matrix scaled by 5 (multiplied by 5) so tiiat it contains integer 
valued elements which simplifies multiplying by * e^^ in the LMS update equation. 



W 



k+1 = ^ 



Wk-H-PviV2^-(X,.eK)- 



(21) 



The scaling factor is taken into account in the selection of the gain value \l. Constraining the gain to g and tiie phase to 
kjc at the normalized frequency of 0.5 simplifies implementing the orthogonal projection matrix PviV2-'-: half of the ele- 
ments are zero and the otiier half are either +1,-1, or +4. Thus, multiplying the projection matrix Pv^ V2-^ by tiie gradient 
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\^ues X|(*^ requires only shift registers and adders. 

The ACQ/TRK signal shown in FIG. 8B disables adaptation of the FIR filter during acquisition, that is. while acquir- 
ing the acquisition preamble 68 shown in FIG. 2B. Thus, the adaptive equalizer B103 adapts only after acquiring the 
acquisition prean^e 68. 

5 

Reduced Cost Orthogonal Constraint Matrix 

Even though the above orthogonal projection matrix Pv^ V2-'- has a simple structure wherein half of the elements are 
zero, rt may not be cost effective to directly inplement it due to the significant number of shift and accumulate operations 
10 necessary to compute PviV2''"*(X,^*e,^).ln order to reduce the cost and connplexity, an alternative embodiment of 
the present invention decimates the modified LMS adaptation algorithm C1 1 4 as illustrated in FIG.s 1 1 A. 11 B, 11 C and 
11D. 

A mathematical basis for the circuit of FIG. 1 1 A will be discussed before describing the details of operation. Refer- 
ring again to the above equation (20), 
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Pv/ = |.Pv, = l-V,(v/v,)-\/ 



by combining the above and V2 vectors into a Nx2 matrix 
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(those skilled in the art will recognize that shifting the time base provides four alternatives for the V matrix) then the 
30 operation (V/VJ'^ of equation (20) reduces to 



35 



1 0' 
51 
0 5 



I' 



Thus, equation (20) reduces to 



Pv-^ = l-1/5-W^. 



40 



Multiplying both sides of equation (22) by 5 provides 

5»Pv-^ = 5-l- W. 
46 Referring again to equation (21), 

WK.i = WK-n-Pviv/.{X,.eK) 
setting X ,^ • e ^ = g ,^ reduces equation (21) to 



50 



WK-M=Wk-^-Pv^-9K = W,-n.[5-g,-w^-gJ. 



Defining 6 = v • g ,^ 



(22) 



(23) 



(24) 
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r+1 0-10 ...1 

[o -1 0 +1 ...J 



§0 

83 



(25) 



w 



IS 



20 



25 



that is. 5o=go-g2+ -- and 6 ,=-9 1+93- 
then computing • 



+1 0" 






+5o 




0 -1 












-1 0 






-So 


(26) 


0 +1 




h 




+5; 





and computing 5 • g ,^ - • g provides 



30 


+5o 




5go-5o 




-5, 














35 


+5, 







{21) 



40 

FIG. 1 1 A implements the above equations (25), (26), (27) and urtimately equation (21) in order to update the coef- 
ficients of the equalizer filter ClOO according to the modified LMS adaptation algorithm of the present invention. To fur- 
ther reduce the implementation cost, the circuit of FIG. 1 1 A decimates the adaptation algorithm by the number of filter 
coefficients (10 in the example shown); that is, the adaptation algorithm operates only on every tenth sample value 
45 rather than on every sanple value and updates only one f liter coefficient per clock period. This is illustrated by expand- 
ing equation (24) 



50 



55 



w 




W 

^ (k + 0)0 




5g(kH. 0)0-5 


W 

^ (k + 2) I 




^(k + l)l 




5g(k + l)l+5 


^(k + 3)2 




^(k + 2)2 




5S(k + 2) 2'*"5 


W 

^(k+4)3 




^(kH.3)3 




5g(k + 3) 3-5 



(28) 



where k+l Is the clock period and 9 (qj = • e j . Thus, at each clock period, the next gradient value g(i)j can be com- 
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puted by multiplying a sample value latched every tenth clcx:k cycle by the next error value ej. The new gradient g^j is 
then used to update the corresponding filter coefficient VJ^-^ of equation (28). 

Referring now to the circuit shown in FIG. 1 1 A. the sample values are input into the discrete time equalizer filter 
CI 00 and the equalized sample values input into the interpolator B122 simitar to FIG. 8B. A serial-tOi)arallel cirojit 

5 C140 converts the interpolated sanrple values Yk+^ into even and odd subsequences Y2k+x and Y(2k.i)+T. where the 
notation 2k indicates that two new interpolated sample values are output from the serial-to-parallel circuit CI 40 at every 
other sample period. A slicer C142 generates the conresponding even and odd estimated subsequences ~Y2k+T artd 
*Y(2k.i)+T which are subtracted from the interpolated sample values at respective adders (C144^.C1442) to form an 
even and odd sample error sequences e2k.K and B^^k-iy^x- ^ arror value interpolation circuit C146, similar to that of 

10 FIG. 8B, generates the even and odd sample en-or sequences e2k and e(2k-i) which are synchronized to the A/D sanrple 
rate. 

As mentioned above, the circuit of FIG. 11 A decimates the adaptation algorithm of the present invention by the 
number of coefficients in the equalizer filter CI 00. For the example shown, the equalizer filter C100 has 10 filter coeffi- 
cients; accordingly, a decimate by 10 circuit CI 48 loads a sample value into a delay register CI 50 every tenth sample 

75 period. Thereafter, the output of delay register CI 50 is represented by X2k-2j-i where: j=1 ->5 incremented by 2 at every 
other sannple period. The output of the delay register CI 50 X2k-2j-i 'S multiplied by the sample en-ors e2k and e(2k-i) at 
respective multipliers (CI 52^. CI 522) to form tiie gradient values g2j and g2j.i used in equation (24). 

The gradient values g2j and g2j-i are then shifted into respective shift registers (CI 54i .CI 542). To implement equa- 
tion (27), the gradient values g2j and g2j.i are multiplied by respective alternating ±1 (C156i,C1562) and accumulated 

20 in respective accumulators (C158i.C1582). After accumulating 5 gradient values in each accumulator (C1 58^, CI 582). 
the outputs of tiie accumulators (which represent 60 and of equation (25)) are latched by respective decimate by 5 
circuits (C160i,Cl602) and the accumulators (C158i.Cl582) are reset. The values Sq and 61 are then multiplied by 
respective alternating ±1 (C162i,C1622) to implement equation (26). The gradient values g2j and g2j.i at the outputs of 
the shift register (Cl54i,C1542) are multiplied by 5 (C164i.Cl642) and the values 60 and 5i are subtracted tiierefrom 

25 at adders (C166-I.C1662) in order to implement equation (27). 

To finish tiie adaptive update algorithm (i.e., to implement equation (28)), tiie output of adders (CI 66.), CI 662) are 
scaled by a gain factor \x (C168^.C1682) which is reduced by a factor of 5 to account for the scaled up projection oper- 
ator. The output of tfie gain factor n (C168i,C1682) is subtracted at adders (C170i.C1702) from tiie con-esponding filter 
coefficient (W2j.i,W2^ selected by a multiplexor (C1 72^ .CI 722) from a bank of registers (C174i,Cl742). The ADAPT2j 

30 signal selects the 2p coefficient from tiie bank of registers (C174^,C1742) for updating. After subtracting the update 
value, the updated filter coefficient (C176i,C1762) is restored to the bank of registers (C174i,C1742) and used by tiie 
equalizer filter C100 during the next clock period to equalize the sample values according to its updated spectrum. 

Decimating the update algorithm by the number of filter coefficients, as in FIG. 1 1 A, decreases the implementation 
complexity but at the cost of slowing convergence of the equalizer filter CI 00 toward an optimum setting. This is 

35 because the decimated update algoritiim does not use all of the sample values to compute the gradients gj. FlGs. 1 1 B 
and 1 1C show an alternative embodiment of the decimated update algorithm of the present invention wherein more of 
the sample values are used to compute the gradients gj. which adds complexity but improves performance because tiie 
equalizer filter C100 converges faster. 

The circuit in FIG. 1 1 B operates similar to the circuit in FIG. 1 1 A described above except for the addition of respec- 

40 tive gradient averaging circuits (Cl76i.C1762) which compute an averaged gradient g2j.i and g2j over several sample 
values. 



45 n = 0 



where N is a predetermined number of sample values out of the number of equalizer filter taps. In one embodiment, all 
of the sample values could be used (i.e. , a number of sample values equal to the number filter taps), or in an alternative 

50 embodiment, a decimated number of sample values could be used in order to reduce the complexity and cost. 

FIG. 1 1C shows an embodiment of the gradient averaging circuits (C176^,C1762) wherein tiie number of sample 
values is decimated by 2; that is, 5 out of the 10 sample values corresponding to the 10 equalizer filter taps are used to 
compute the averaged gradients g2j.i and g2j. 

In operation, a decimate by 2 circuit CI 78 of FIG. 1 1 B stores every other sample value in the delay register CI 50. 

55 The output of the delay register CI 50 is multiplied (CI8O1-I8O5) by the en'or value e2k and delayed versions of tiie en-or 
value e2k-n» and tiie results are accumulated in accumulators (C182.,-1825). After accumulating five gradients, the con- 
tents of tiie accumulators (C182^-1825) are transferred to registers (C184^-184s), and the accumulators (C182^-1825) 
are cleared. Then, at every other sample period, ttie contents of registers (C1 84^-1845) are shifted from left to right (i.e., 
CI 845 = CI844; CI844 = CI843; etc.) and tiie output of register CI 845 is the averaged gradient g2j output by the gra- 
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10 



dient averaging circuit (0176^.01762). The averaged gradients Q21 and ggj.! (CI 86^.01862) are then used to update the 
coefficients of the equalizer fitter CI 00 using circuitry CI 88 in the same manner as described with reference to FIG. 
11 A. 

FIG. 1 1 D illustrates yet another embodiment of the present invention which further reduces the cost and conplex- 
ity, as compared to FIG. 11 A. by updating the even and odd filter coefficients sequentially. That is, during the first N peri- 
ods (where N is the number of filter coefficients) the circuit of FIG. 11 D computes the coefficient update values of 
equation (28) for the even filter coefficients (Wq, W2, W4,-.), and then during the next N sample periods it computes the 
update values for the odd filter coefficients (W^, W3. W5....). A decimation circuit CI 90 decimates the error value e^^ by 
two, and approximately half the circuitry as that of FIG. 1 1 A is used to compute the update values. TTie decimate by two 
circuit CI 90 is actually syncopated; that is, it outputs the error values for k=0, k=2, k=4, k=6, k=8, and then outputs the 
error values for U^, k=1 1. k=13. k=15, k=17 (assuming the adaptive filter CI 00 comprises ten filter taps). 

IVIathematically, operation of the update circuit in FIG. 1 1 D can be described by. 
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assuming the adaptive filter CI 00 comprises 10 filter taps. 

Obviously, the embodiment of FIG. 11 D will decrease the performance of the adaptive filter C100 due to the 
decrease in convergence speed. However, the gradient averaging circuit of FIG. 110 can be used to improve the per- 
formance of the circuit in FIG. 11 D, similar to the embodiment of FIG. 1 1 B. Thus in the embodiments of FIG. 1 1 A-1 1 D, 
the performance versus cost and complexity varies-the preferred configuration is selected according to the require- 
ments of the user. 

Although the interpolated timing recovery and adaptive equalizer of the present invention have been disclosed in 
relation to a d=0 PR4 read channel, the principles disclosed herein are equally applicable to other types of sampled 
amplitude read channels including d=1 EPR4 or EEPR4 read channels. In a d=1 read channel, for example, the slicer 
B141 of FIG. B4A is replaced by a pulse detector as described in the above reference U.S. Patent No. 5,359,631. 

Furthermore, the particular constraint frequency of 1/4Ts used In the disclosed embodiment is not intended to be 
limiting. Other constraint frequencies could be used without departing from the scope of the present Invention. For 
example, a 3T preamble could be used In which case the constraint frequency would be 1/6Ts (if constraining to the 
preamble frequency). 

Still further, those skilled In the art will appreciate the many obvious modifications that are possible to the adaptive 
equations disclosed herein. For example, a shift in the time base would change the FIR filter's magnitude and phase 
response at the constraint frequency which would result in different V-t and V2 matrices. Also, the modified LMS algo- 
rithm, 

WK.i = WK-Ji-PviV2'--(XK-e,) 
could be implemented after rearranging terms, 

Wk*i = WK-H-eK-(Pv,V2^.X,) 



19 



EP0 805 447 A2 



or 

Wk.i=WK-M-XK-(Pv,V2^-e,). 

These, arxJ like modrfcations. are within the scope of the present invention. 

The objects of the invention have been fully realized through the embodiments disclosed herein. Those skilled in 
the art will appreciate that the various aspects of tiie invention can be achieved through different embodiments without 
departing from the essential function. The particular embodiments disclosed are Illustrative and not meant to limit the 
scope of the invention as appropriately construed by the following dainr^. 



Table 1 



Channel 


Transfer Function 


Dipulse Response 


PR4 


(1-D)(1+D) 


0,1,0, -1,0, 0, 0. ... 


EPR4 


(1-D)(1+D)2 


0. 1. 1. -1. -1.0. 0. ... 


EEPR4 


(1-D)(1+D)3 


0, 1.2. 0, -2. -1.0,... 
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Table B2 



5 



W 



15 



20 



25 



35 



45 



50 



T*32/Ts 


C(-2) 


C(-2) 


C(0) 


c(1) 


C(2) 


C(3) 


0 


0.0000 


-0.0000 


1.0000 


0.0000 


-0.0000 


0.0000 


1 


0.0090 


-0.0231 


0.9965 


0.0337 


-0.0120 


0-0068 


2 


0.0176 


-0.0445 


0.9901 


0.0690 


-0.0241 


0.0135 


3 


0.0258 


-0-0641 


0.9808 


0.1058 


-0.0364 


0.0202 


4 


0.0335 


-0.0819 


0.9686 


0.1438 


-0.0487 


0-0268 


5 


0.0407 


-0.0979 


0.9536 


0.1829 


-0.0608 


0-0331 


6 


0.0473 


■0.1120 


0.9359 


0.2230 


-0-0728 


0.0393 


7 


0.0533 


-0.1243 


0.9155 


0.2638 


-0-0844 


0.0451 


8 


0.0587 


-0.1348 


0.8926 


0.3052 


-0.0957 


0.0506 


9 


0.0634 


-0.1434 


0.8674 


0.3471 


-0-1063 


0.0556 


10 


0.0674 


-0.1503 


0.8398 


0.3891 


-0.1164 


0.0603 


11 


0.0707 


-0.1555 


0.8101 


0.4311 


-0.1257 


0.0644 


12 


0.0732 


-0.1589 


0.7784 


0.4730 


-0-1341 


0.0680 


13 


0.0751 


-0.1608 


0.7448 


0.5145 


-0-1415 


0.0710 


14 


0.0761 


-0.1611 


0.7096 


0-5554 


-0-1480 


0.0734 


15 


0.0765 


-0.1598 


0.6728 


0.5956 


-0.1532 


0.0751 


16 


0.0761 


-0.1572 


0.6348 


0-6348 


-0.1572 


0.0761 


17 


0.0751 


-0.1532 


0.5956 


0-6728 


-0.1598 


0.0765 


18 


0.0734 


-0.1480 


0.5554 


0-7096 


-0.1611 


0.0761 


19 


0.0710 


-0.1415 


0.5145 


0-7448 


-0.1608 


0.0751 


20 


0.0680 


-0.1341 


0.4730 


0-7784 


-0.1589 


0.0732 


21 


0.0644 


-0.1257 


0.4311 


0-8101 


-0.1555 


0.0707 


22 


0.0603 


-0.1164 


0.3891 


0-8398 


-0.1503 


0.0674 


23 


0.0556 


-0.1063 


0.3471 


0-8674 


-0.1434 


0.0634 


24 


0.0506 


-0.0957 


0.3052 


0.8926 


-0.1348 


0.0587 


25 


0.0451 


-0.0844 


0.2638 


0.9155 


-0.1243 


0.0533 


26 


0.0393 


-0.0728 


0.2230 


0.9359 


-0.1120 


0.0473 


27 


0.0331 


-0.0608 


0.1829 


0.9536 


-0.0979 


0.0407 


28 


0.0268 


-0.0487 


0.1438 


0.9686 


-0.0819 


0.0335 


29 


0.0202 


-0.0364 


0.1058 


0.9808 


-0.0641 


0.0258 


30 


0.0135 


-0.0241 


0.0690 


0.9901 


-0.0446 


0.0176 


31 


0.0068 


-0.0120 


0.0337 


0-9965 


-0.0231 


0.0090 



55 Claims 

1. A sampled amplitude read channel for reading digital data from a sequence of discrete time sample values gener- 
ated by sampling an analog read signal from a read head positioned over a magnetic medium, comprising: 



21 



EP0 805 447 A2 

(a) a sampling device for sampling the analog read signal to generate the discrete time sample values: 

(b) an adaptive equalizer, responsive to the discrete time sannple values, for generating equalized sanople val- 
ues according to a target response: 

(c) interpolated timing recovery for generating interpolated sample values: and 

5 (d) a discrete time sequence detector for detecting the digital data from the interpolated sample values. 

2. The sanripled amplitude read channel as recited in claim 1 . wherein: 

(a) the adaptive equalizer comprises an FIR fitter cx>mprising a plurality of filter coefficients: 

10 (b) the adaptive equalizer operates according to a least mean square algorithm responsive to an en'or value e^^ 

computed as a function of an output of the FIR filter and an estimated ideal value: 

(c) the adaptive equalizer conrtprises an interpolation circuit, responsive to the interpolated sample values, for 
generating the en^or value e^ synchronous with the discrete time sample values; and 

(d) the least mean square algoritiim uses the en^or value e^^ to adjust the filter coefficients of the FIR filter. 

15 

3. Hie sampled amplitude read channel as recited in claim 2, wherein the error value e^ is computed relative to a dif- 
ference between an output of tiie interpolated timing recovery and a sample value estimator responsive to the inter- 
polated sample values. 

20 4. The sampled amplitude read channel as recited in claim 2. wherein the error value e^ is computed relative to a dif- 
ference between an output of tiie FIR filter and an output of the interpolation circuit. 

5. The sampled amplitude read channel as recited in claim 2, wherein: 
25 (a) the least mean square algorithm is: 

(b) Wk are tiie filter coefficients of the FIR filter; 
30 (c) n is a predetermined gain; and 

(d) is a discrete time sample value. 

6. TTie sampled amplitude read channel as recited in claim 2. wherein the adaptive equalizer further comprises a con- 
straint circuit for constraining a frequency response of the adaptive equalizer at a predetermined constraint fre- 

35 quency. 

7. The sampled amplitude read channel as recited in claim 6. wherein tiie constraint circuit comprises an orthogonal 
projection operation for constraining the filter coefficients relative to the predetermined constraint frequency. 

40 8. The sampled anplitude read channel as recited in claim 7, wherein: 

(a) the constraint circuit operates according to a modified least mean square algoritiim, 

Wk.i = WK-n-PviV2-^.(X,.e,) 

45 

(b) Wk are tiie filter coefficients; 

(c) ^ is a predetermined gain; 

(d) Py/iV2^ is tiie orthogonal projection matrix; and 

(e) Xk is a discrete time sample value. 

50 

9. The sampled amplitude read channel as recited in claim 1 , wherein the interpolated timing recovery further gener- 
ates a data clock for clocking the discrete time detector. 

10. T?ie sampled amplitude read channel as recited in claim 1. wherein: 

55 

(a) the analog read signal comprises pulses modulated by the digital data at a predetermined baud rate; and 

(b) the sampling device samples tiie analog read signal at a rate above tiie baud rate. 

11. The sampled amplitude read channel as recited in claim 1. wherein the interpolated timing recovery comprises an 
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interpolator responsive to an interpolation interval t proportional to a time difference between a channel sannple 
value and a desired Interpolated sample value. 

The sampled amplitude read channel as recited in claim 1, wherein the analog read signal comprises pulses mod- 
ulated by the digital data at a predetermined baud rate, the interpolated timing recovery further comprising: 

(a) a frequency offset generator for generating a frequency offset Af proportional to a frequency difference 
between a sampling dock and the baud rate; and 

(b) a mod-Ts accumulator for accumulating, modulo-Ts. the frequency offset Af to generate an Interpolation 
interval t where Is is a predetermined sample period of the sampling clock. 

The sampled anplitude read diannel as recited in claim 12, wherein the frequency offset generator comprises: 

(a) a phase error detector for detecting a phase en^or AO between an interpolated sample value and an esti- 
mated sanple value; and 

(b) a loop filter for filtering the phase error Ae to generate the frequency offset Af. 

The sampled amplitude read channel as recited in claim 11. wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k) which approximates an ideal impulse response of an ideal inter- 
polation filter: 

sinc(7t • (k-x/Ts)); where: 

k is a time index; and 
25 Ts is a sample period of the sampling clock. 

15. The sampled amplitude read channel as recited in claim 11, wherein the interpolator connprises an actual interpo- 
lation filter having an actual impulse response h(k) generated by minimizing a mean squared error between a fre- 
quency response of the actual interpolation filter and an ideal frequency response corresponding to an ideal 

30 interpolation filter. 

16. The sampled amplitude read channel as recited in claim 1 1 , wherein: 

(a) the interpolator conprises an actual interpolation filter; and 

35 (b) a plurality of coefficients ^ *he actual interpolation filter are computed in real time as a function of x. 

17. A method for reading digital data from a sequence of discrete time sample values generated by sampling and ana- 
log read signal from a read head positioned over a magnetic medium, comprising: 

40 (a) sampling the analog read signal to generate the discrete time sample values; 

(b) adaptively equalizing the discrete time sample values to generate equalized sample values according to a 
target response; 

(c) interpolating the equalized sample values to generate interpolated sample values; and 

(d) detecting the digital data from the interpolated sanple values. 

45 

18. The method for reading digital data as recited in daim 17, wherein the step of interpolating sample values com- 
prises the step of generating a data dock for clocking the discrete time sequence detector. 
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